Bosonic Mode Mixing in the Superconducting State Spectral Function of 
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Photoemission spectra of Bi2Sr2CaCu208+i below Tc show two features near the (7r,0) point of the zone: 
a sharp peak at low energy and a higher binding energy hump. We find that the sharp peak persists at low 
energy even as one moves towards (0, 0), while the broad hump shows significant dispersion which correlates 
well with the normal state dispersion. We argue that these features are naturally explained by the mixing 
of electrons with a bosonic mode which appears only below Tc, and speculate that the latter may be related 
to the resonance seen in recent neutron data. 
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Angle resolved photoemission data on the quasi-two 
dimensional high temperature superconductors can be 
interpreted in terms of the one-electron spectral func- 
tion [y. This implies that important information about 
the self-energy E, and how it changes from the normal 
to the superconducting (SC) state, can, in principle, be 
obtained by analysis of the ARPES lineshape. This obvi- 
ously has important ramifications in elucidating a proper 
microscopic theory of high temperature superconductors. 

Perhaps the most dramatic effect in this regard is the 
temperature dependence of the lineshape in Bi2212. A 
very broad normal state spectrum near the (tt, 0) point of 
the zone evolves quite rapidly for T < Tc into a sharp, res- 
olution limited, quasiparticle peak ||l[ followed at higher 
binding energies by a dip [0J3] then a hump, the latter 
corresponding to where the spectrum recovers to its nor- 
mal state value. Similar effects are observed in tunneling 

i- 

In this paper we focus on another remarkable difference 
between the normal state and SC state data which has 
not been noticed earlier. In Fig. 1, we show spectra for 
a Tc =87K Bi2212 sample along T-M-Z, i.e., (0, 0) - 
(7r,0) — (27r,0), (a) in the normal state (105 K) and (b) 
in the SC state (13 K) from which we note two striking 
features. First, we see that the low energy peak in the SC 
state persists over a surprisingly large range in k-space, 
even when the normal state spectra have dispersed far 
from the Fermi energy. For example, the sharp peak is 
visible at about 40meV even in curve 4 of Fig. lb, when 
the corresponding normal state spectrum is peaked 320 
meV below Ep. Second, when the hump in the SC state 
disperses, it essentially follows that of the normal state 
spectrum. 



To highlight these two features, we plot in Fig. 2 the 
position of the sharp peak and that of the hump in the SC 
state versus the single (broad) spectral peak position in 
the normal state. The experimental points in this figure 
are, by themselves, highly suggestive of the quantum me- 
chanical mixing of two levels in the SC state electronic 
spectral function. We note the similarity of this plot 
to the spectrum of an electron interacting with a sharp 
phonon mode ||5|] in conventional superconductors. At the 
same time it is very important to stress that, in marked 
contrast to the phonon case, the "bosonic mode" seen 
here appears only in the superconducting state, and, as 
we shall argue below, has a many-body origin. 

We will discuss below a phenomenological self-energy 
for the SC state which incorporates the mixing of electron 
states with a bosonic mode. The resulting spectral func- 
tions are not only able to capture the features of the line- 
shape for a given k, but, more importantly, also able to 
describe the very different dispersions of the sharp peak 
and the broad hump. Note that, while near kp the sharp 
peak is quasiparticle-like in nature, it is mostly bosonic- 
like once the higher energy hump has begun to disperse. 
It is likely that this bosonic mode has the same physical 
origin as that recently observed by neutron scattering in 
YBCO [|. 

The data of Fig. 1 were obtained on high quality 
slightly overdoped Bi2212 single crystals {Tc = 87 K), 
with measurements carried out at the Synchrotron Radi- 
ation Center, Wisconsin, using a high resolution 4m nor- 
mal incidence monochromator |Q. 22eV photons polar- 
ized along T — M (the Cu-0 bond direction) were used for 
both narrow energy scans (resolution FWHM^IS meV) 
and wide energy scans (FWIIM=35 meV). Similar results 



were seen on a variety of samples with different doping 
levels, photon polarizations, and photon energies. 

We have already noted above the persistence of the low 
frequency peak in the SC state from curves 4 through 
16 in Fig. lb, about half of the Brillouin zone. As one 
moves from either extreme towards M (curve 11), (1) the 
broad hump at high binding energy moves towards Ep 
and (2) weight is transferred from the hump to the low 
frequency peak, which is is fairly fixed in energy. The 
same phenomena are also seen along M to Y (Fig. Ic). 

A potential complication in interpreting these data is 
the influence of images of the CuO bands (ghosts) caused 
by the incommensurate superlattice ||3|,[7|, since one pre- 
dicts a Fermi crossing of one of these images near curve 
4. The following arguments can be made against such 
an effect. First, the ghosts are not visible in the normal 
state in this polarization geometry. They do however be- 
come quite visible if the photon polarization is rotated by 
45° 0. Moreover, comparison of superconducting state 
spectra in these two polarizations indicate that the mid- 
point of the leading edge in the present polarization (20 
meV) is near that of the M point, whereas in the 45 de- 
gree rotated polarization, the midpoint is 5 meV. The 
latter value is consistent with the expected gap size on 
the ghost band for this k value. Second, the intensity at 
the peak position (35 meV) monotonically rises from F 
with a maximum near Af , indicating only one spectral 
feature. 

We can also eliminate the two bilayer-split bands as the 
explanation of the non-trivial lineshape with a peak, dip 
and hump. First, we recall our earlier results on the 
intensities of the peak and the hump at M as a function 
of varying photon polarization direction from in to out of 
plane. The fact that the intensities of both these features 
scaled together argues against these arising from two dif- 
ferent bands. Second, the very different dispersions of 
the peak and hump noted in this paper also implies the 
same: if the sharp peak was a second band it must also 
show up above Tc, but it does not. 

We begin our discussion of self-energy effects by re- 
calling earlier work on the SC state lineshape. At low 
energies, there is a linewidth collapse [g[ upon cooling 
through Tc, leading to resolution limited spectral peaks 
[Q at T <C Tc. This is most easily understood (at least 
near kp) by the freezing out of electron scattering once 
the SC gap opens up. That is, the electron lines in the 
Feynman diagram shown in Fig. 3a are gapped. At higher 
energies (3A for the s-wave case) the imaginary part of 
S (JttiE), which determines the linewidth, must go back 
to its normal state value, thus leading to a qualitative 
explanation M of the dip and hump in SC state spectra 
near kp. For a gap with nodes, relevant to the cuprates, 
these arguments have to be modified, since contributions 
to IttiT, exist all the way to zero frequenc y pO[ . 

More recently, theoretical calculations JL1[, motivated 
by the neutron scattering experiments in YBCO , have 



found a sharp structure in the SC state particle-hole sus- 
ceptibility (bubble in Fig. 3a). Basically, the dominance 
of Q = (tt, tt) scattering, coupled with the effect of the 
sign change of the d-wave gap with Q on the BCS coher- 
ence factors of the bubble, leads to a sharp resonance at 
ujo between A^ax and 2Amax [Amax = A(^_o))- In some 
models there is a true collective mode at frequency loq- 

Motivated by the data in Fig. 2 and the above consid- 
erations, we incorporate the effect of a bosonic mode on 
the SC state spectral function, taking proper account of 
both the real and the imaginary parts of the self energy. 
We note that the importance of structure in ReYl has not 
been fully appreciated in the context of high Tc super- 
conductors. A mode at loq leads to a a singularity in the 
SC state Iml] at uiq + A. By Kramers-Kronig transfor- 
mation, this implies a singularity in /JeS at loq -\- A. This 
feature in /JeS in turn can lead to an additional pole in 
the spectral function. A simple way to visualize this is to 
to think of an electron interacting with an Einstein mode 
at Wo . One obtains for this model the well-known dis- 
persion with two branches resulting from the nrixing of 
the electron and the mode, with in general one pole be- 
ing primarily electron-like in nature and the other being 
primarily mode-like. In addition, there is no damping for 
frequencies lower than ujq + A. 

We model the SC state data using the phenomenolog- 
ical self energy Il3] 

/mS(tj) = ToN{\lo\) + TiNQiol ~ wq), |w| > wq + A 

= roiV(|w|), A<|w|<wo + A 

= 0, |a;|<A (1) 



where N{iX}) = Lojyjuj'^ — A^ is the BCS density of states, 
and 



7ri?eE(w) = ro7V(-tj) In | - ^ + \J uP- - A'^\l A 



-hriiV(wo - w) In |wo - w -H ^J(uj - wo)2 - A2|/A 

-{w^-w} (2) 
The spectral function is then [^ 



7rA(w) = Ira 



Z2(w2-A2)-e2 



(3) 



with (a complex) Ziuj) = 1 — Yiiuj^jio. 

Note that at large w the model has a constant linewidth 
broadening, so that we are effectively ignoring the lo- 
dependence of the normal state JmE. We have found 
that simulations with a large constant IrnS, which also 
take into account background emission, can adequately 
fit the normal state data (see also [131). From nor- 
mal state data, we find that Fi is of order 200 meV 
along T — M. The constant linewidth simplification al- 
lows us to directly obtain the dispersion et from tight 
binding fits to the normal state peak position M. Our 



final results are not very sensitive to Fq which is in- 
cluded for numerical stability; we use 30 meV, but sim- 
ilar results are found for 5 meV. (To obtain signifi- 
cant damping of the low frequency pole would require 
making A complex.) For A, we assume a d-wave gap 

Ak = A„a^(cOs(fca;a) - C0s{kya))/2 with A,„aa; = 32 

meV. The best agreement with experiment is found by 
choosing the bosonic mode frequency loq — l.SAk- 

The resulting real and imaginary parts of S at M are 
shown in Fig. 3b emphasizing the singular behaviors at 
A due to the Fq term and at wq + A due to the Fi 
term. In Fig. 4, we show a comparison of the result- 
ing spectral function (convolved with the experimental 
energy and momentum resolution) to experimental data 
at M for both wide and narrow energy scans. To aid in 
the comparison to experiment, a step edge background 
(which gives a good model of data for k > kp) is added 
to the calculation. The positions of the sharp peak and 
the hump obtained from the calculations are compared 
to the experimental data in Fig. 2. The relative lack of 
movement of the low frequency peak is well reproduced, 
as well as its lack of visibility for curves 1 through 3 in 
Fig. lb. This lack of movement is due to several factors: 
(1) the fiatness of the band dispersion Ck near M, (2) 
the lowering of loq ^ Ak as one moves towards F, and 
(3) the infiuence of both ReY, and the gap. The last is a 
new effect worth commenting on. The self-energy (Eq. 2) 
implies a mass enhancement [Z > 1) in the SC state rel- 
ative to the normal state, which acts to push spectral 
weight towards Ep. On the other hand, A itself pushes 
spectral weight away from Ep. Thus the dispersion is 
dramatically fiattened relative to the normal state. 

We now discuss the p dependence of E, which was es- 
sentially ignored, except in so far as it was crudely cap- 
tured by that of wq- Assuming that Q = (tt, tt) scattering 
dominates, as suggested by the neutron data, we would 
replace FiiV(a; -|- ujq) in Eq. 1 by g^p^gAp+Qiu + uq) 
(for a; < 0) where g is the interaction vertex. Assuming 
a quasiparticle pole approximation for A when solving 
Eqs. 1 and 2, this would imply a dip in the spectrum 
at |a;| = Ep+g + ujq where E^ = el/{ReZp)'^ + A^, and, 
as before, a persistent low frequency peak if Z is large 
enough. Note that p and p + Q are coupled in the equa- 
tion for E which requires a numerical solution. Although 
we leave this complexity for future work, we note that 
this coupling implies that if a low frequency peak exists 
for p, then one also exists for p + Q (in fact, they self- 
consistently generate one another if g is large enough). 
This is just the effect observed in the data along (tt, 0)- 
(tt, tt), shown in Fig. Ic, in that a persistent low fre- 
quency peak exists for about the same angular range as 
that along (tt, 0)-(0, 0). Thus the expected Q dependence 
of the scattering is indeed refiected in the ARPES data. 
We note that such scattering can also cause a pairing 
instability in the d-wave channel. 

In conclusion, we have shown the presence of a per- 



sistent low frequency peak in photoemission spectra in 
Bi2212 in the SC state which exists over a large angular 
range near the M point. The dispersion of this feature 
and the higher binding energy hump as a function of mo- 
mentum indicates that the electrons in the SC state are 
interacting with a bosonic mode of resonant character 
with a frequency between A^ax and 2A,inax- Our results 
once again emphasize that the self-energy is dominated 
by electron-electron interactions, which is consistent with 
an electron-electron origin to the pairing. 
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FIG. 1. EDCs in (a) the normal state (105 K) and (b) the 
superconducting state (13 K) along the line V — M — Z, and 
(c) the superconducting state (13 K) along the line M — Y , for 
a slightly overdoped {Tc — 87K) Bi2212 sample with photon 
polarization F — M. The zone is shown as an inset in (c) with 
the curved line representing the observed Fermi surface. 



FIG. 2. Positions (eV) of the sharp peak and the broad 
hump in the SC state versus normal state peak position. Solid 
points connected by a dashed line are the experimental data, 
the solid lines are obtained from the calculations described 
in the text, and the dotted line represents the normal state 
dispersion. 



FIG. 3. (a) Feynman diagram for the lowest order con- 
tribution to the self-energy from electron-electron scattering. 
The solids lines are electrons, the dashed lines the interaction, 
(b) JmE and ReZ at M from Eqs. 1 and 2 for the parameter 
set given in the text. 



FIG. 4. Comparison of the data at M for (a) wide and 
(b) narrow energy scans with calculations based on Eqs. 1-3, 
with an added step edge background contribution. 
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Figure 2 Norman et al. 
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Figure 3 Norman et al. 
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Figure 4 Norman et al. 



